Quantum chaos in A =46— 50 atomic nuclei 
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In recent years many authors have shown great interest in the fluctuation properties 
of energy levels. It is well known that the fluctuation properties of quantum systems with 
underlying classical chaotic behaviour and time-reversal symmetry agree with the pre- 
dictions of the Gaussian Orthogonal Ensemble (GOE) of random matrix theory, whereas 
quantum analogs of classically integrable systems display the characteristics of the Poisson 
statistics [1]. 

One of the best systems for the study of quantum chaos is the atomic nucleus, which 
was studied in pioneer work [2] and has been the subject of many investigations [3]. In 
atomic nuclei, the fluctuation properties of experimental energy levels are best studied in 
the domain of neutron and proton resonances near the nucleon emission threshold, where a 
large number of levels of the same spin and parity in the same nucleus are present, and an 
excellent agreement with GOE predictions has been found [2]. In the ground state region, 
however, the samples of consecutive experimental levels of the same spin and parity in any 
one nucleus are quite small. Therefore it is more difficult to calculate reliable mean values 
and fluctuations of statistics such as energy level spacings. In order to circumvent this 
difficulty, in recent years [4-6] statistical analyses of experimental low energy levels have 
combined data from a large range of excitation energies and angular momenta of a nucleus 
or a set of nuclei. Such analyses have provided evidence suggesting that spherical nuclei 
show level spectra close to GOE predictions and deformed nuclei show strong deviations 
from GOE behaviour. In a recent analysis [7] of level spacings close to the yrast hne 
of deformed nuclei with Z — 62-75 and A — 155-185, the average level spacing for 
states with the same spin and parity was calculated for the total ensemble instead of for 
individual nuclei. The level spacing fluctuations obtained are quite close to the Poisson 
distribution, showing evidence of regular motion, although it is possible that this Poisson 
behaviour may be related to a hidden symmetry in the K quantum number. In general 
it seems quite clear that the available experimental data are insufficient to establish the 
borderlines in mass number, excitation energy, etc., between order and chaos in nuclei. 
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Theoretical work, especially realistic shell model calculations, should help to estabhsh 
the domains of chaos in nuclei. The statistical analyses of shell-model energy spectra 
and wave functions have been mainly concentrated on the sd shell region, and a very 
chaotic behaviour has been found for these nuclei [8-11]. However, recent shell-model 
calculations by Bae et al. [12] have shown that for nuclei of mass A — 212 shell-model 
spectra display features of regular motion in some cases. The main reason for this seems 
to be the relatively small values of the residual interaction matrix elements as compared 
to the average spacing between neighboring single-particle levels in heavy nuclei. 

In this work, we undertake the statistical analysis of the shell-model energy levels in 
the A — 46-50 region. Exact calculations are performed in the (/7/2,P3/2,/5/2,Pi/2) shell- 
model space, assuming a ^'^Ca inert core [13]. The diagonalizations are performed in the 
m-scheme using a fast implementation of the Lanczos algorithm with the code ANTOINE 
[14] . For a fixed number of valence protons and neutrons we calculate the energy spectrum 
for projected total angular momentum J and total isospin T. The interaction we use is a 
minimally modified Kuo-Brown realistic force with monopole improvements [15]. 

We calculate the T — states from J = to J = 9 for all the combinations of 6 
active nucleons, i.e. ^^V, ^^^Ti, ^^Sc, and ^^Ca, and also for ^^Ca and ^^Ca. Table 1 shows 
the m-scheme dimensions and the largest dimensions in the JT scheme of the isotopes 
analyzed. 

Since we are looking for deviations from chaotic features, we are mainly interested in 
the low-lying levels, up to a few MeV above the JT yrast fine. Let us consider the energy 
levels up to 4, 5 and 6 MeV above the yrast line, and calculate the fiuctuations around the 
average spacing between neighboring levels. In this range of energies, the level spectrum 
can be mapped into unfolded levels with quasi-uniform level density by using the constant 
temperature formula. In order to guarantee that the results up to different energies are 
unaffected by the unfolding procedure, the unfolding is performed using always the whole 
set of levels up to 6 MeV, for each JT set in the nucleus. In some cases, the resulting 
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number of levels is too small, e.g. 3 levels for J = in ^^Ca, and then we use a larger set 
of levels for the calculation of the mean level spacing as a function of energy. The mean 
level density can be assumed to be of the form 

-p{E) = ^e^v[{E-E,)/Tl (1) 

where T and Eq are constants. For fitting purposes it is better to use not p{E) but its 
integral N{E). We write 

N{E) = r p{E')dE' + TVo = exp [{E - E^)/T\ - exp [-E^/T] + TVq- (2) 
Jo 

The constant A^o represents the number of levels with energies less than zero. Following 
Shriner et al. [5], we consider Eq. 2 as an empirical function to fit the data and let A^o 
take non-zero values. The parameters T, £"0 and A^o that best fit N{E) are obtained by 
minimizing the function: 

G{T, Eo, No) = r^'^^[N{E) - N{E)fdE, (3) 

J E 

where N{E) is the number of levels with energies less than or equal to E. The energies 
E^in and E^ax are taken as the first and last energies of the level sequence. As an 
example, Fig. 1 illustrates the fit to the integrated level density N{E) for the J'^T = 6+1 
levels of ^^Ti. 

The spectral statistic P{s) is used to study the local fluctuations of the energy levels 
[16,17]. P{s) is the distribution of nearest-neighbour spacings Sj = N{Ei+i) — N{Ei) of 
the unfolded levels. It is obtained by accumulating the number of spacings that lie within 
the bin (s, s + As) and then normalizing P{s) to unity. 

For quantum systems whose classical analogs are integrable, P{s) is expected to follow 
the Poisson hmit, i.e. P{s) — exp(— s). On the other hand, quantal analogs of chaotic 
systems exhibit the spectral properties of GOE with P{s) — (7r/2)s exp (— [1,2]. 

The distribution P{s) is the best spectral statistic to analyze shorter series of energy 
levels and the intermediate regions between order and chaos. In order to quantify the 
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chaoticity of P{s) in terms of a parameter, it can be compared for example to the Brody 
or the Berry-Robnik distributions, which are adequate for the description of intermediate 
situations between order and chaos. Although each of these distributions have some 
advantage in limiting cases [18], they are very similar in a particular case like ours. We 
use here the Brody distribution [19], given by 

P{s, uj) = a{(jj + 1)5"^ exp (4) 

with 

a^{r[^]r\ (5) 

This distribution interpolates between the Poisson distribution {cu — 0) of integrable 
systems and the GOE distribution {cu — 1) of chaotic ones, and thus the parameter cu can 
be used as a simple quantitative measure of the degree of chaoticity. 

In order to obtain more meaningful statistics, P{s) is calculated using the unfolded 
level spacings of the whole set of J = 0-9 levels for fixed T up to a given energy limit 
above the yrast hne. The number of J = 0-9 spacings below 4, 5 and 6 MeV range from 
42, 66 and 105 in ^^Ca, to 86, 149 and 231 in ^^Ti, respectively. 

Table 2 shows the best fit Brody parameter uj of the P{s) distribution for the J = 0- 
9 set of level spacings in the A = 46 nuclei up to 4, 5 and 6 MeV above the yrast 
line. Clearly, ^^V, ^^Ti and ^^Sc are chaotic for these low energy levels, but there is 
a considerable deviation from GOE predictions in ^^Ca, which is a single closed-shell 
nucleus. In view of the peculiarity of this nucleus, we performed calculations for ^^Ca 
and ^°Ca, and obtained again strong deviations toward regularity, as the values of Table 
2 show. Thus, for the Ca isotopes we find the same kind of phenomenon obtained by Bae 
et al. [12] in the heavy single-closed nuclei ^^^Rn and ^^^Pb, namely that low-lying states 
deviate strongly from chaoticity toward regularity. 

To obtain a better estimate of the Brody parameter, we can combine spacings of 
different nuclei. In Fig. 2 (a) we plot P(s) for ^^V+^^Ti+^^Sc, and in Fig. 2 (b) for 
^^CaH-''^Ca+^°Ca, up to 6 MeV above the yrast lines. The number of level spacings is 



now sufficiently large to yield meaningful statistics and we see that Ca isotopes are not 
very chaotic at low energy, in contrast to other nuclei in the same region. 

Concerning the energy dependence of uo in the low energy region, for individual nuclei 
we see fluctuations which may be due to the fact that some of the samples of levels are 
not very large. Table 2 shows the results for ^^V+^^Ti+^^Sc and ^^Ca+^^Ca+^^Ca. For 
these larger sets of levels, fluctuations are reduced, and the energy dependence of uj is 
seen to be small in the 4-6 MeV range, although the chaoticity increases slightly with the 
excitation energy. 

Since the semi-closed nuclei are quite regular for the low-lying levels, compared with 
other neighboring nuclei, it is interesting to find out in the framework of the model the 
chaoticity at higher excitation energies. Therefore we have calculated the whole spectrum 
for ^^Ca, unfolding the levels with the local density method [3]. Different J values were 
combined to improve the P{s) statistics. The value of uj increases rather smoothly with 
the excitation energy end point, and we obtain uj — 0.72, 0.83 and 0.88 for levels below 
12, 16 and 20 MeV, respectively. We do not find significant differences in uj for different 
J values in ^^Ca when they are analyzed separately. 

Why are Ca isotopes less chaotic than their neighbors? We observe that the two-body 
matrix elements of the proton-neutron interaction are, on average, larger than those of 
the proton-proton and neutron-neutron interactions. Consequently the single-particle 
mean-field motion in nuclei with both protons and neutrons in the valence orbits suffers 
more disturbance and is thus more chaotic. 

It should be noted that an analysis [20] of experimental energy levels below 4.3 MeV 
excitation energy in the semi-magic nucleus ^^^Sn yields a near-neighbor spacing distri- 
bution which is intermediate between GOE and Poisson, with uj — 0.51 ±0.19. This result 
is consistent with the theoretical findings of Bae et al. [12] for ^^^Rn and ^^^Pb, and our 
present results for Ca isotopes. 

The energy spectra of all these nuclei are the result of the interaction of like nucleons. 
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and we conclude that this residual interaction is relatively weak and, at least for low 
energy levels, it only partially destroys the regular motion of nucleons generated by the 
nuclear mean field. 

We can ask ourselves why do all shell-model calculations give chaotic features for 
sd shell nuclei, without any significant deviations towards regularity. First, it should be 
noted that most of these calculations include a large number of states, up to excitation 
energies far above the nucleon emission threshold. This is, for example, the case of the 
^^O calculations of Bae et al. [12], which include the full set of levels for several J values 
and obtain lu = 0.96. Second, we notice that Ormand and Broglia [10] obtained a GOE- 
like distribution for the first two spacings of each spectrum for a set of sd shell nuclei. 
However, these nuclei have valence protons and neutrons and are thus similar to the 
case of ^^V, ^^Ti and ^^Sc, for which we also find chaotic behaviour. We conclude that 
no regular features have been found in the sd shell region because single-closed nuclei 
have a very small configuration space for these nuclei and thus too few low-lying levels 
for statistical analysis, and also because the disturbance of single-particle motion by the 
two-body interaction is greatest in light nuclei, where the range of the single-particle 
orbits is not much larger that the range of the nuclear force. 

This work has been partially supported by MURST and INFN, Italy, and by DGICYT 
grant PB93-0263 and CICYT grant AEN-0776, Spain. 
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Table 1 

Matrix dimensions in the m scheme and maximum dimensions in the JT scheme for the 
analyzed isotopes. 





46Y 


46Ti 


^^Sc 


^^Ca 


48Ca 


50Ca 


m-scheme dimension 


121 440 


86 810 


30 042 


3 952 


12 002 


17 276 


largest dim. J+T 


4+0 


4+1 


4+2 


4+3 


4+4 


4+5 


J~^T dimension 


4 750 


8 026 


3 783 


615 


1 755 


2 468 
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Table 2 



Brody parameter uj for the nearest neighbour level spacings distribution for < J < 9, 
T — Tz states up to 4, 5 and 6 MeV above the yrast hne in the analyzed nuclei. 



Energy 


46v 






^6Ca 


^^Ca 


50Ca 


46v+46rpi^46Sc 


46Ca+48Ca+50Ca 


< 4 MeV 


1.14 


0.90 


0.81 


0.41 


0.58 


0.67 


0.92 


0.56 


< 5 MeV 


1.10 


0.81 


0.96 


0.53 


0.58 


0.69 


0.93 


0.60 


< 6 MeV 


0.93 


0.94 


0.99 


0.51 


0.66 


0.62 


0.95 


0.61 
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Figure Captions 



Figure 1: Best fit of tlie constant temperature formula to the integrated level density for 
J^T = 6+1 of ^^Ti. 

Figure 2: P{s) distribution for low-lying levels of fp shell nuclei with < J < 9: (a) 
^^V, ^*5Ti and ^^Sc; (b) ^^Ca, ^^Ca and ^"Ca. The dotted, dashed and sohd curves stand 
for GOE, Poisson, and Brody distributions, respectively. 
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